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Abstract. The Polyakov loop is an order parameter to the confinement-deconfinement phase transi- 
tion in pure gluodynamics. Recently a simple phenomenological model has been proposed to descri- 
be the available lattice data for the renormalized Polyakov loop in the deconfinement phase |1]. 
These data exhibit unequivocal inverse power temperature corrections driven by a dimension two 
gluon condensate. The heavy quark free energy is also studied and compared to lattice data. 

Introduction. In pure gluodynamic the confinement-deconfinement transition can be 
characterized by the breaking of the Z{Nc) discrete symmetry, belonging to the center of 
the gauge group, SU(A^c)- The order parameter is the traced Polyakov loop, defined by 

L{T) = (tr a{x)) = (^\r P (^e'^ lo'^ dxoAo{x,xo)'^ ^ ^ 

where ( ) denotes vacuum expectation value. Aq is the gluon field in the (Euclidean) time 
direction. A perturbative evaluation of the Polyakov loop was carried out in [2]. 

The Polyakov loop correlation functions are related to the change in the free energy 
arising from the presence of a static quark-antiquark pair in a thermal medium. The 
colour singlet free energy yields the heavy quark potential at finite temperature, 

Fi(x, T)/T = -log(TrL(x)Lt(0)) +c{T) . (2) 

A renormalization procedure for the Polyakov loop based on the computation of 
singlet and octet correlation functions between Polyakov loops in the limit of large 
separation has been considered in recent lattice studies II^IU]. The additive normalization 
constant, c{T), is ambiguous and can be fixed at short distances by comparison with the 
zero temperature heavy quark potential, i.e. for the singlet free energy: Fi{r -C 1/7") — 
Vqq{f)- With such a prescription the renormalized Polyakov loop is defined as yllj]: 

In Ref. [1], we have proposed a model to describe the available lattice data for U™. Here 
we will show that it also describes consistently the lattice results for the free energy. 
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Phenomenological model for the Polyakov loop. From now on we work in the 
Polyakov gauge: (9oAo(x, jcq) = 0. We assume that, in the deconfinement phase, the field 
Ao(x) is sufficiently well described by a Gaussian distribution. From ([T]) one finds ^ 



L = exp 



(4) 



To describe the dynamics of the Ao(x) field we use the 3-dimensional reduced effec- 
tive theory of QCD Let Doo(k)5ai, denote the 3-dimensional propagator, then 

{aL)t = {nI-\)T j ^^oo(k) . (5) 

In perturbation theory the propagator becomes D^q{^) = 1/ (k^ + t^D) ("^o ~ T is the 
Debye mass), and when inserted in (O it reproduces the known perturbative result 
to LO We will take into account the non-perturbative contributions coming from 
dimension -2 condensates. So, we consider adding to the propagator new phenomenolo- 
gical pieces driven by positive mass dimension parameters: ^ 

<(k) = 7r^^- (6) 
Adding the perturbative and non perturbative contributions, we get for the Polyakov loop 

Comparison to lattice data. The perturbative contribution to the Polyakov loop has 
logarithmic corrections in T. We will assume that (Aq ^)^^ is constant in T, or at most 
logarithmic. This means that logL has a temperature power correction. 

The lattice data of the renormalized Polyakov loop of Ref. [3] for Nf = flavors and 
Ref. 10] for Nf = 2 flavors have been displayed in Fig.[T] These data suggest a linear fit 

of the form -21ogL = a + b{Tc/Tf which yields 

0.23(1), A^/ = 2 2 .yvp_/(0.87±0.02GeV)2, A^^ = 



-0.26(1), Nf = 2' 6 Y^o,a/i |(0.71±0.01GeV)2, Nf = 2 

with x^/DOF = 0.45, 6.14, for Nf = 0,2 respectively. The perturbative results for the 
highest temperature 6Tc are in qualitative agreement with the fitted values for a, Eq. ^ 
Finite temperature results for the pressure in pure gluodynamic ^ yield for the 
gluon condensate (0.93(7) GeV)^ (in the temperature region used in our fits and in 
Landau gauge). Our result is also in reasonable agreement with existing zero temperature 
values ilO]. 



^ The Gaussian ansatz is correct up to &{g^) in perturbation theory, and exact in the large Nc limit. 
^ Such ansatz parallels those make at zero temperature in the presence of condensates 0]. 
We consider the LO and NLO perturbative result of Ref. |3l 
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FIGURE 1. The logarithmic dependence of the renormalized Polyakov loop versus the inverse tempera- 
ture squared in units of the critical temperature. Lattice data from |3, 4]. At the left we plot fits with 
a adjustable constant or predicted by NLO perturbation theory [2]. Purely perturbative LO and NLO 
results for A'/ = are shown for comparison. The right figure illustrates the duality between the Polyakov 
loop and the quark-antiquark potential at zero temperature. The line represents rV^g(r), takes from lattice 
data L121 . and modified with the change r = l/nio- 



Heavy quark free energy. The quark-antiquark potential can be related to the scat- 
tering amplitude corresponding to one gluon exchange. In the non-relativistic limit 

F,{x,T)^-^-g'l^^e'^-^Doo{l^). (9) 

We can consider non perturbative contributions for the free energy by adding the non 
perturbative term ® to the perturbative propagator, i.e. £)oo(k) = ^ooC*) +^0' C^)- 
LO, ff{g^), and NLO, &{g^), the singlet free energy has the form 

' ' INc \A%rN^-\ T ) 2Nc ^ IN J ' 

(10) 

Taking the value obtained in ^ for g^{A^ 0)7^^ match the free energy to lattice 

data in]. This way, one obtains the r and T dependence of the running coupling a^. 
This is displayed in Fig. Q. One can see a smooth behaviour, and it confirms the better 
fit of the Polyakov loop with a = constant. The values are relatively small, which is in 
contrast with existing analysis of at finite T lisLITlll. 

In we consider the limit of large separation in the free energy, Eq. dlOb . one reobtains 
Foc(r) ~ —2T\ogL. In the zero temperature limit (we take nior ^ in (fTOb ) 



\t2 1 2 a3/42 \NP 

^7r: + — r = ^m{r). (H) 



which is just the quark-antiquark potential at T = 111211 . The Coulomb term is the 
standard perturbative result at LO. The second one is a non perturbative linear contri- 
bution, and we get for the string tension o = g^iA^ ^j^q/2Nc- From Eqs. ^ and (flTT) 
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FIGURE 2. The r dependence of the running couphng in pure gluodynamic for different values of 
temperature. It is obtained by mathing Eq. ( I10> to lattice data for the free energy from Ref. fllll . 



we deduce a property: 

F^iT)=Vqq{r)l^,/^^, (12) 

which is valid if we assume that a.,(r,r = 0)(A^^)^£o = a^(r °°,7^)(^o,fl)r^- This 
duality between F^{T) and Vqq{r) is only valid at LO in perturbation theory. ^ 

To check numerically Eq. Jl^ . we plot in Fig.[l]the lattice data for —IblogL versus 
can see a remarkable agreement. This duality looks deeply into the analogy 
between the quark-antiquark potential at zero temperature and the Polyakov loop. 
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^ The duality is formal, i.e. we are assuming g = constant ^ g{r, T). If we take into account the different 
assymptotic behaviours of a^: o:j(r) = as(r, T = 0) and (r) = as{r °o,T); Eq. il2\ writtes 

bF^iT)=Vqc,{r)\r = r/T, with = (a,(r)/a,(r))3/4 and 7= (a,(r)/a,(r)3)i/^ (13) 



